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Abstract. Riemannian first-passage percolation (FPP) is a continuum model, with a distance function 
arising from a random Riemannian metric in M d . Our main result is a shape theorem for this model, which 
says that large balls under this metric converge to a deterministic shape under rescaling. As a consequence, 
we show that smooth random Riemannian metrics are geodesically complete with probability one. 



o 

1. Introduction 

1.1. Motivation. We introduce Riemannian first-passage percolation (FPP) as a model of random geometry 
in the continuum. Our work is motivated by standard FPP, where one builds a random distance function on 
Z d from independent, identically-distributed "passage times" over bonds of the lattice. We instead construct 
' our distance function on W 1 from a random Riemannian metric. The main result of the paper is a shape 

p I . theorem for Riemannian FPP: if Bt is the random Riemannian ball of radius t, then \B t tends toward 

a limiting shape. From this it follows that a random Riemannian metric is geodesically complete with 
probability one. 

Hammersley and Welsh [7] introduced standard FPP in 1965 in order to model fluid flow through porous 
media. Consider the d-dimensional lattice 1 d with d > 2. Let {tb} be a family of independent, identically 
distributed, non-negative random variables, indexed by bonds (nearest-neighbor edges) b of the lattice. For 

. any z, z' £ Z rf , define 

>; = inf£>, 

(X; , 7 f-^ 

be-y 

where the infimum is taken over all lattice paths 7 connecting z to z'. This d is a random distance function 
on Z d . For a very good introduction to standard FPP, see Howard [8] . 
[ — ■ Consider the random distance d(0, nei) between the origin and the point nei = (n, 0, . . . , 0). One wishes 

to study the asymptotic behavior of this quantity as n — > 00. In |12j . Kingman formulated his famous 
subadditive ergodic theorem in order to prove the basic result of FPP: provided the passage times have finite 
mean, there exists a non-random constant fi ei , such that 



lim id(0,nei) = ^ c 



almost surely and in L 1 . The same is clearly true for all coordinate axes and, more generally, for each 



direction v £ S d 1 , there exists a non-random constant [i v , such that 
(1.1) lim ±-d(0,nv) = fx v 

almost surely and in L 1 , where nv £ 7L d is the nearest lattice point to nv. The constant /i„ is non-zero 
provided that the probability that tb = is less than the critical percolation probability for 1 d [5]. 
The shape theorem of Cox and Durrett [2] is a stronger result. Consider 

B t = {z £ Z d : d(0,z) < t}, 

the random ball of radius t in Z d . This is a lattice object, so we "inflate" it to get a continuum one: for 
z £ Z d , let C z = [z — 1/2, z + l/2) d be the unit cube centered at z in R d , and let 

B t = (J C z . 

zeB t 

We define the rescaling jBt as the set of all points x £ M. d such that tx £ Bt- The shape theorem says 
that there exists a non-random convex, compact set A, depending only on the distribution of tb, such that 
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jB t — > A: for all e > 0, with probability one, there exists a time T such that if t > T, then 

(l-e)AC \B t C (l + e)A 

Just as standard FPP is a model of random geometry on the discrete lattice Z d , Riemannian FPP is a 
model of random geometry in the continuum. We consider a random Riemannian metric g on M. d whose 
distribution is translation-invariant, has finite-range dependence and satisfies certain moment conditions. 
By the standard construction in Riemannian geometry, this defines a random distance function d(x, y) in 
M. d . Kingman's theorem can again be applied to prove that for each unit vector v, there exists a non-random 
constant /j, v > such that ^d(0,nv) — » fj, v a.s. and in L 1 . From the positive-definiteness of g, we prove 
Theorem 12.51 [i v > for all v £ S^ -1 . Consider the set 

A = {xeR d :\x\<n-} M }, 

and the Riemannian ball of radius t centered at the origin 

B t = {x£R d : d(0,x) < t}. 

Theorem 13. II is the shape theorem for Riemannian FPP, which states that for all e > 0, with probability one, 
there exists a random time T > such that if t > T, then 

(l-e)AC \B t C (l + e)A 

Consequently, we call A the limiting shape of the model. 

Theorem [313] follows from the shape theorem: with probability one, no curve 7, parametrized by Riemann- 
ian length, reaches infinity in finite time. When the metric is further assumed to be smooth with probability 
one, then this is geometrically significant: by the Hopf-Rinow theorem |13) . this is equivalent to geodesic 
completeness of the metric. 

To prove our results, we need some technical estimates on the distance function d, which we obtain in 
Section [21 To prove these, we discretize the continuum model: to each point z £ Z d , we associate a certain 
value X z based on the Riemannian metric g(x) over the unit cube C z = [z - 1/2, z+ l/2) d . By treating X z 
as defining a dependent FPP model on the lattice, we prove some estimates of the "entropy-energy type" on 
X z , and from these derive the desired estimates on d. 

1.2. Geometry Background and Notation. Before introducing any probabilistic structure, we introduce 
some geometric notation. Consider R d with d > 2 and the standard Euclidean coordinates. Write 

SPD = {symmetric, positive-definite d x d real matrices}, 

and let g £ C(M. d , SPD) be a continuous matrix-valued function on W 1 with values in SPD. g defines a 
Riemannian structure on R d : for tangent vectors v, v' £ T x M. d , we consider the inner product (v, g(x)v'}. For 
a single vector v, we denote by ||w|| = y/ (v, g{x)v) and \v\ = yf (v, v) the Riemannian and Euclidean lengths 
of v, respectively. For a C 1 -curve 7 : [a, b] — > R d , we define the Riemannian and Euclidean lengths of 7 by 

R{i) = I Il7(*)ll di and L( 7 ) = / |<y(i)| dt, 

respectively. We say that a curve is finite if it has finite Euclidean length; for our model, Theorem 13.41 will 
imply that finite curves have finite Riemannian length. The Riemannian distance between two points x and 
y is defined by 

d(x,y) = inf i?( 7 ), 

7 

where the infimum is over all C^-curves 7 connecting x to y. 

For a Riemannian metric g, we define the real, positive functions 

A(x) = maximum eigenvalue of g{x) and M x ) = minimum eigenvalue of g{x). 

For any K C R d , define 

A(K) = sup A(x) and X(K) = inf X(x). 

By the continuity and positivity of g, if K is bounded then 

< \{K) < A(K) < x). 



A SHAPE THEOREM FOR RIEMANNIAN FIRST-PASSAGE PERCOLATION 



3 



For z € Z d , let C z = [z — 1/2, z + l/2) d be the unit cube centered at z. Write 

A z = A(C.) and \ z = \{C Z ). 

1.3. Riemannian FPP. Let il = C(M. d , SPD) and let T be the cr-algebra generated by cylinder sets. Let 
P be a translation-invariant probability measure on which has finite-range dependence, and consider 

a random Riemannian metric g G Q with distribution P. The finite-range dependence means there exists 
some R > such that if \x — y\ > R, then g{x) and g(y) are independent. Furthermore, suppose that Ao has 
a finite moment-generating function. That is, 

(1.2) M(r) = E[e rAo ] < oo for all rel, 

where E denotes expectation with respect to P. By translation invariance, the family {A 2 } is identically 
distributed, and we refer to its generic element as A; similarly for {X z } and A. By Chebyshev's inequality 
[4], (|1.2p implies that A and A have exponential tail decay: 

P(A > u) < P(A > u) < M(r)e _r ", 

for all u > and r > 0. 

We provide a concrete example. 

Example 1.1. Let c : [0, oo) ->lbea compactly-supported covariance function (see [5] for examples). Let 
£ : R d -> K be a mean-zero, stationary, istropic Gaussian field with covariance function c; that is, 

E[£(z)] = and E[Z(x)t{y)] = c(\x - y\) 

for all W d . The covariance c must satisfy certain necessary and sufficient conditions [17] for the field 

£ to be everywhere continuous with probability one; suppose this is the case. 
Let g : M. d — > K be the diagonal matrix with entries 

ff«(aO=log(l + e«*>), 

for 1 < i < d. This is continuous and positive, so it suffices to show that the assumption ()1.2j) is satisfied. 



Proposition 1.2. Assumption (jl.2l) is satisfied for this choice of g, and the shape theorem (Theorem 13. If) 
applies. Let A be the limiting shape, defined in (13.11) . The measure P is isotropic, so Corollary 13.21 implies 
that A is a Euclidean ball. Furthermore, if the field £ is C 1 with probability one, then Corollary 13.51 implies 
that the metric g is geodesically complete. 

Proof. Assume for simplicity that c(0) = 1. The Gaussian concentration inequality [1] implies that 
i > u) = P(sup<e > log(e u - 1)) < cxp(- log(c u - l) 2 /2) < 2e-" 2 / 2 . 



The fundamental theorem of calculus and Fubini's theorem imply that 

(/■Ao \ roo poo 

l + J re ru duj=l + J rc™P(A > u) du < 1 + 2 J re ru e~ u2 /2 du, 

which is finite for all r. Thus (II. 2j) is satisfied, and the results of this paper apply to the random Riemannian 
metric g. □ 

The random variables X z and A z give rise to a dependent FPP model on sites of the lattice Z d . In Section 
12.11 we prove some general estimates for dependent FPP, then in Section l2~2l we apply these to estimates on 
our distance function d. Our techniques are based on the energy-entropy methods of mathematical physics, 
where one shows that an event occurs with extremely low probability over one particular connected set 
("high energy"), but sums this over all possible connected sets at the origin ("high entropy"). One adjusts 
parameters in the problem so that this sum converges, then applies the Borel-Cantelli lemma. A large- 
deviations estimate like 11.21 is critical: the number of connected sets at the origin grows exponentially in n, 
the size of the sets, so the probabilities must decay exponentially in n for the arguments to hold. 

In the standard FPP setting of passage times tb across bonds 6, Cox and Durrett [5] prove that a necessary 
and sufficient condition for a shape theorem is that Emin{£i, . . . ,t2d} d < oo, where fcj are 2d independent 
copies of Thus we believe that our assumption (|1.2j) is not the most general, and can be replaced by a 
finite moment estimate on A instead to prove a more general result. 
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In probability theory, Kingman's subadditive ergodic theorem [3] is used to prove that stationary, subad- 
ditive sequences obey laws of large numbers. If X n ^ m is a non-negative, stationary sequence which satisfies 
X n .m < X nir +X ri7n , then ^^o,n converges almost surely and in L 1 . Furthermore, if the sequence is ergodic, 
this convergence is to a non-random constant. In our context, the sequence in question is X njn = d(nv, mv) 
for a fixed unit vector v. The subadditivity condition is exactly the triangle inequality for d, so Kingman's 
theorem implies that for each v £ S x , there exists a non-random constant fi v > such that 

(1.3) lim ±d(0,tv) = fi v 

almost surely and in L . The constants /i„ may depend on the direction v, though if the measure P is 
isotropic (rotationally-invariant) then = /i will be independent of v. We show in Theorem 12 . 5 1 that /i„ > 
for all v. 

Proposition 1.3. fi v is a continuous function of v. 

Proof. This remarkably short proof is due to Kesten [TT] (see his Proof of Theorem 1.7 on page 158). First 
we show that the function fi v is bounded above as a function of 3 d " 1 . Write v = Yli v%e ii f° r the standard 
basis vectors e^ in R d . We use the triangle inequality to bound d(0, tv) by the sum of the distances between 
successive points 0, tv e%, tv 1 e\ + tv 2 e2 and so on until tv. Translation-invariance implies 

Ed(0, tv) < Ed(0, tw x ei) + • • • + Ed(0, tv d e d ). 

Dividing by t and taking the limit t — > oo gives 

d 

fj, v = lim ±Ed(0,tv) < V lim ±Ed(0, tv l e d ). 

i=l 

The terms which equal zero we may ignore; for the non-zero terms, we make the substitution t' — so 
that the right-hand side equals 

d d 
V lim ^Ed(0,±t' ei ) = V l^l/ie, < dmax{^ c J, 

— ' t'— >oo ^ — ' 

i=l i=l 

as desired. 

Now, consider two different unit vectors v and v', and write u = njE^jn ■ By the same arguments, 

\[i v — < lim jKd(tv,tv') = (i u \v — v'\. 
This tends to zero as v' — > v since fj, u is bounded above. □ 

2. Discretization Lemmas 

2.1. Dependent FPP on a Lattice. For a continuous curve 7, we would like to introduce a discrete 
analogue T on the lattice. For example, z G T if 7 meets the cube C z . However, this set T is not connected 
on Z d ; consider the straight line from to (1,1) in R 2 . We get around this by modifying the familiar 
graph structure of Z d to introduce a new lattice, which we call the ^-lattice. In this section, we prove some 
estimates for dependent FPP on the *-lattice, then in Section l2~2l we apply these estimates to the continuum 
model. 

For z £ Z d , we write z — (z 1 , . . . , z d ). We say that z, z' £ X d are ^-adjacent if maxi<i< £ ;(z — z') 1 < 1. 
The *-lattice is the graph with vertex set Z d , and edge set given by ^-adjacency; that is, the usual lattice 
Z d along with all the diagonal edges. 

We say that a set T C 7L d is ^-connected if for all z, z' £ T, there is a path from z to z' along the ^-lattice 
which remains in the set T. Technically, that there is a finite sequence of ^-adjacent points beginning with 
z and ending with z' , all contained in T. 

Let S n be the number of ^-connected sets which contain the origin. 

Lemma 2.1. There exists a such that S n < er™. Obviously, a > 1. 
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Proof. Clearly, log S n is a non-negative subadditive sequence: 

log S n+m < log S„ + log S m . 

By Fekete's lemma [5], there exists a constant a > such that log S n < an for all n. Defining a = e a proves 
the result. □ 

Let X z be a stationary, non-negative random field on the ^-lattice with finite-range dependence, and with 
a finite moment-generating function 

(2.1) M(r) = E[c rX ] < oo for all rel. 

The finite-range dependence means that there is an integer R > 1 such that if \z — z'\ > R, then X z and X z < 
are independent. 

If r C Z d is a collection of lattice points, we write 

x(r) = J> z , 

and call this the passage time of T. 

The following two lemmas can be thought of as spatial laws of large numbers. The first says that for 
sufficiently large n, if there is a uniform bound on X(T)/n, then there is also a uniform bound on |T|/n. 
The second lemma reverses the implication, though with different constants. 

Lemma 2.2. Suppose X z additionally satisfies 

(2.2) F(X = 0) < a-^ R+1 ^ d . 

For any A > there is a non-random B > such that, with probability one, for any sequence a n € Z d , there 
exists N > such that for all n > N, if V is a ^-connected set which contains the point a n and X(T) < An, 
then |r| < Bn. 

Lemma 2.3. For any B > there is a non-random C > such that, with probability one, for any sequence 
a n E Z d , there exists N > such that for all n > N, if T is a ^-connected set which contains the point a n 
and |r| < Bn, then X(T) < Cn. 

Some assumption like (|2.2I) is necessary for Lemma \2. 2 1 Let p = ¥(X — 0), and suppose that p > p c , the 
critical probability for site percolation on the ^-lattice. By percolation theory, with probability one, the set 

T = {z eZ d : X z = 0} 

contains an infinite ^-connected component. That is, |T| = oo but X(T) = 0. No such assumption is 
necessary for Lemma l2.3l 

In this paper, our assumptions of non-negativity and (|2.2j) are stronger than necessary: we apply these 
lemmas only to the positive fields A 2 and A 2 . However, we anticipate these lemmas to be of independent 
use in future work on Ricmannian FPP, where one may consider fields X z which take the value on a cube 
C z with small but non-zero probability. For example, if E z is the event that A z < h for a sufficiently large 
value of h, one may apply these lemmas to X z = 1e,, the indicator function of E z . 

The generality of the sequence a n is needed for Lemma 12.41 In the proof of that result, we fix a point 
x G M. d , and define the sequence a n — nx £ Z d to be the nearest lattice point to nx. 

Proof of Lemma \2.2\ In what follows, we assume that T is a ^-connected set. Fix A > and the sequence 
a n . Consider the events 

E n = {3 T such that a n € T, X(T) < An, and |T| > Bn}. 

We claim that we can choose an integer B large enough so that ¥(E n ) decays exponentially in n. From the 
Borel-Cantelli lemma it will follow that, with probability one, only finitely many of the events E n occur, 
which will prove the result. 

Let Zi be a predetermined enumeration of Z d ; for example, a spiral path beginning at a n . For any "it- 
connected set r, we define f C T by proceeding along the sequence Zi and including each point of T at a 
distance at least R away from the previous points chosen. In the form of an algorithm: 
• Let ii be the first index for which z,, e T. Let gi = z^. 
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• Given {fll,.. . jgj— 1}, let ij be the first index for which Zj. £ T and so that \zi j — gji\ > R for 
1 < f < 3- Let gj = z ir 
Let r' = {31,52, • ■ • }• If r is finite, then so is V . Let 

B(z, R) = {z' &Z d :\z- z'\ < R} 

be the Euclidean ball of radius R centered at z in Z d . Note that \B(z, R)\ < (2R + l) d . To avoid this 
cumbersome factor (2R + l) d which appears frequently, for the remainder of this proof we write 

K = (2R+l) d . 

By construction, the set T is covered by taking balls around every point in L': 

LC |J B(z,R). 

Let r be a ^-connected set described in the event E n , so that 

Bn<\T\ < \B(z,R)\ < \T'\K. 
zer' 

Thus, |r'| > Bn/K. 

We may assume that T consists only of the first Bn points it meets of the sequence Zi ; the non- negativity 
of X z implies that the passage time X(T) is still at most An. Similarly, we assume |r'| = \_Bn/K\, the 
integer part of Bn/K. Furthermore, since fCT, 

X(T') < X(T) < An. 

Thus the probability ¥(E n ) is bounded by 

(2.3) P (3 F s.t. a n G F, X(T') < An, and |r| = Bn) <J2 F ( X ( T ') ^ An ) > 

r 

where the outer sum is taken over all ^-connected sets T for which a n € T and |F| = Bn. Lemma 12 . 1 1 implies 
that the number of such sets is bounded by a Bn . 

The family of random variables {X z }r' is independent, since the points z € V are separated by distances 
at least R. Let {Xi} be \Bn/K~\ independent copies of X. The exponential Chebyshev inequality [3] implies 
that the right-hand side of (|2.3|) is bounded above by 

/[Bn/K] \ 

a Bn p Xi < An < a Bn e rAn E (e~ r ^ Xi ) 

for any r > 0. Again by independence, if we write M(—r) = Ec~ rX , this is equal to 

a Bn e rAn( Ee -rXyBn/K] ^ ^ n M( _ r )fln/ K = f &T A/B M( _ r) l/K^ S " _ 

By the bounded convergence theorem, as r tends to infinity, M(—r) — > P(X = 0), which is strictly less 
than a~ K by assumption (|2.2I) . Let r be large enough so that M(—r) < a~ K . Write p = <jM(—r) x l K < 1, 
so that 

HEn) < ( P C rA ' B ) Bn . 

Choose the number B to satisfy 

(2.4) rAj log(l/p) < B < rAj log((l + p)/2p). 
This implies 

i±2 < pe rA / B < 1. 

The left inequality will be used later in the proof of Thereom 12.51 the right inequality implies that 

< $>e^/ fl )*™ < oo, 

so by the Borel-Cantelli lemma, with probability one, only finitely many of the events E n hold. □ 
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Proof of Lemma \2.3\ In what follows, we assume that T is a ^-connected set. Fix an integer B > and the 
sequence a n . Consider the events 

E n = {3T such that a n e T, X(T) > Cn, and |r| < Bn}. 

We claim that we can choose C large enough so that F(E n ) decays exponentially in n. From the Borel- 
Cantelli lemma, with probability one, it will follow that only finitely many of the events E n occur, which 
will prove the result. 

As in the proof of Lemma |2~21 we consider passage times of subsets of T to exploit independence. Unlike 
in that proof, it does not suffice to consider just one subset T' . Instead, we partition the set into k disjoint 
subsets r 1; . . . , Tfe, and consider the passage times of each. 

There are k — R d points in the cube {0, . . . , R — l} d ; order them z\, . . . , Zk- We can partition the lattice 
Z d into k subsets by considering i?-translations of these points. We write z = (z 1 , . . . , z d ) for all z € Z d . For 
any *-connected set T, we partition it into k subsets by defining 

Tj = {z € T : R divides z % — for alH = 1, . . . , d} 

for j — 1, . . . , k. Since points in Tj are separated by distance at least R, for each j the family of random 
variables {A z }r is independent. 

Let r be a set described in the event E n . We may assume that |T| = Bn, since if it is less, the inclusion 
of additional points will only increase the passage time calculation. Since 

X(T)=X(T 1 ) + ---+X{r k ), 

if X(T) > Cn, then X(Tj) > Cn/k for some j. Furthermore, < |T| < Bn for each j. 
As in Lemma \2.2l the probability V(E n ) is bounded above by 

(2.5) 5^p(a j e{l,...,fc} s -t- x ( t j) > Cn/k), 

r 

where again the sum is over ^-connected sets V for which a n € T and |T| = Bn. For each j, the family 
of random variables {Xi}rj is independent. Let {X{\ be Bn independent copies of X. The exponential 
Chebyshev inequality [3] implies that the right-hand side of (|2.5p is bounded above by 

/ Bn \ 

k<j Bn ¥ [J2 Xi > Cn / k ) < ka Bn e- Cn / k E(e^ Xi J , 

where we have added extra independent copies of Xi to make the number of terms exactly Bn. If we write 
M = M(l) = Ec x , then this is equal to 

ka Bn e -Cn/k M Bn = fc ^B e -C/k M B\ n 

We choose O 1 so that a B &~ c l k M B < 1. Thus 

J2 F (E n ) < kY,(v B e- C/k M B Y < oo, 

so by the Borel-Cantelli lemma, with probability one only finitely many of the events E n hold. □ 

2.2. Applications to Continuum Model. In this section, we apply the estimates from the previous section 
to the continuum model. Lemma 12.41 says that at large scales, the Riemannian distance function between 
two points is bounded by a uniform constant K times the Euclidean distance between them. Theorem 12.51 
says that u v is positive for all v. 

The cubes C z form a partition of R d . For x G R , let x be the unique point on the lattice such that 
i£ Cj. The metric g induces two notions of "passage time" over a discrete set T: 

A(7) = ^A, and A(t) = X>. 

Lemma 2.4. There exists a non-random K > such that, with probability one, for all x € M. d and p > 0, 
there exists T(x) > such that if t > T and \x — y\ < p, then 

d(tx,ty) < Ktp. 
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Proof. By rescaling x and y, it suffices to prove the lemma with p = 1. 

Apply Lemma [2.31 with B — 2d and X z = A z . Thus there exists a non-random C > such that, with 
probability one, for any sequence a n E 7L d , there exists N > such that for all n > N, if T is a finite 
♦-connected set which contains the point a n and |T| < 2dn, then A(T) < Cn. 

Let K = C\fd, and suppose that with positive probability, there is some x S M d such that for any n > 0, 
there exists y (depending on n) such that \x — y\ < 1 but d(nx, ny) > Kn. We will show that this leads to 
a contradiction. Let N be as in Lemma 12.31 applied with the sequence a n = nx. 

Suppose that n > N . Let 7 be the straight-line segment from nx to ny. Let 

r = {zez d : 7 nc^0}, 

index the cubes C 2 which 7 meets. Note that nx G T. Clearly, 

|r| < 2d\nx - ny\ < 2dn, 

since \x — y\ < 1. By Lemma l2.31 

A(r) < Cn. 

The distance between nx and ny is a lower bound for the Riemannian length of 7: 

d{nx, ny) < R(j). 

Furthermore, since 7 is a line segment, the Euclidean length of 7 in each cube is at most \fd. Thus we can 
estimate the Riemannian length of 7 by summing A z over F: 

(2.6) d{nx, ny) < i?( 7 ) = ^ R (l n G *) < kz ^ = A ( r )^ ^ Cn ^ = Kn > 

zer zer 

since K = C\fd. This contradicts the assumption that d(nx,ny) > Kn. □ 
Theorem 2.5. The constants fi v are all positive. 

Proof. Suppose fi = fi v = for some unit vector v. 

Let e > 0, and apply Lemma |2~21 with the constant sequence a n = to A = 4e and A z = X z . Thus there 
exists a non-random B > such that, with probability one, there exists iVi > such that for n > JVi, if T 
is a ^-connected set which contains the origin and A(r) < 4en, then |T| < Bn. 

By Kingman's subadditive ergodic theorem, with probability one there exists N 2 > such that if n > N%, 
then 

(2.7) d(0,nv) < en/2, 

since we assumed that /i v = 0. 

Let N — max{iVi, N2}, and suppose n > N. A priori, the distance d(0,nv) need not be realized as the 
Riemannian length of a curve, so let 7 be a C^-curve from to nv with 

R(j) < d(0, nv) + en/2 < en, 

where the second inequality follows from (|2.7|) . 
Define the discrete set 

(2.8) r = {z G Z d : L(jnC z ) > 1/4}; 

that is, z G T provided the Euclidean length of 7 in the cube C z is at least 1/4. 

We claim that T is ^-connected. Suppose not, so that the continuum set W = lj z gr ^» nas a * least two 
components, separated by Euclidean distance at least 1. Let W 1 be the 1/4- neighborhood around W, so 
that the components of W are separated by Euclidean distance at least 1/2. By definition of V, the curve 7 
meets each component of W' , but not the complement R d \ W 1 . Since 7 is continuous, this is a contradiction; 
hence, T is ^-connected. 

In each cube C z , we can estimate the Riemannian length of 7 using A z : 

L^nc z )\ z < Rhnc z ), 

where L denotes Euclidean length. Furthermore, by summing X z over the points of T, we get a lower bound 
for R(<y): 

|A(r) < Lin n c z )\ z < ^ i?( 7 n c z ) < i?( 7 ) < en. 
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Clearly, £ T, so Lemma \2 . 2 1 implies that |T| < Bn. 
In the proof of that lemma, we chose B so that 

B ^ log((l+p)/2p) j4 ' 

for positive constants r and p < 1 not depending on A; see (|2.4p . Since A = 4e, if we write B' = At I log((l + 
p)/2p), then 

(2.9) |r| < B'en. 

Let z be a lattice point in T which minimizes the distance \nv — z\. Clearly, \z\ > n/2. Since T is 
^-connected and z £ T, 

\T\ > n/2\fd. 

For small e, this contradicts (|2.9[) . so /u must be positive. □ 

3. The Shape Theorem and Consequences 

Define the function 

(x) = f/VwM> 

V j [0, x = 0. 

Proposition 11.31 and Theorem 12.51 imply that /i is continuous and, for x ^ 0, strictly positive. It follows from 
the triangle inequality for the distance function that n(x) is a norm on R d . Consider the unit ball in this 
norm, 

(3.1) A = {x : (x(x) < 1} = {x : |z| < At ; / 1 |;i;| }, 
as well as the random Riemannian ball of radius t centered at the origin, 

B t = {x: d(0,x) < t}. 

Theorem 3.1 (Shape Theorem). For all e > 0, with probability one, there exists T such that if t > T, then 

(3.2) (1 - e)A C iS t C (1 + e)A. 

The set A is called the limiting shape of the random Riemannian metric g. The shape theorem for lattice 
first-passage percolation was proved by Cox and Durrett 2 ; our proof is modeled on the arguments in 
Durrett [3]. 

Proof. It suffices to prove the theorem for e £ (0, 1). Let 

(3.3) £ <min {_L i_l±_^ 

where K is as in Lemma \2. 41 Let B E (x,r) denote the Euclidean ball of radius r centered at x. 

We will show that with probability one, there exists T > such that if t > T, then (1 — e)A C \B t . To do 
this, we will first prove that for every x, there is a random T(x) such that if t > T, then the small Euclidean 
ball B E (x, Se 2 ) is contained in \B t . Since the set (1 — e)A is compact, it can be covered by finitely many 
balls B E (xi, 5e 2 ). Letting T — m&x{T(xi)} will prove the result. 

Fix x £ (1 — e)A. We claim that with probability one, there exists T{x) such that if t > T and \x — y\ < 5e 2 , 
then d(0,ty) < t, hence B E (x, Se) C jB t for t > T. By the triangle inequality, 

d{0,ty) < d(0,tx) + d(tx,ty). 

The first term is controlled by Kingman's theorem: with probability one, there exists T\(x) such that if 
t>Ti, then 

d(0, tx) < (1 + e)tp(x) < (1 - e 2 )t, 

since fi(x) < (1 — e). 

The second term is controlled by Lemma T2.4I applied to this x and p = Se 2 . With probability one, there 
exists T2(x) such that for all t > and y with |a; — y\ < Se 2 , then 

d(tx,ty) < KtSe 2 < e 2 t, 

since S < 1/K. 
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Let T(x) = max{Ti,r 2 }. If t > T(x), then for all y £ B E (x,5e 2 ), 

d(0,ty) < (1 - e 2 )t + e 2 t = t, 

implying that B E (x,Se 2 ) C \B t for all t > T(x). Let Xi be finitely many points such that (J B E (xi, Se 2 ) 
covers (1 — e)A, and let T = max^ T(xj). Then for all t >T, 

{l-e)AQ\jB E (x h 5e 2 )clB u 



completing the lower half of the shape theorem. 

Now we prove the upper half. For any x, let T2{x) be as above, so that if t > T2 and \x — y\ < Se 2 , then 
d(tx,ty) < e 2 t. By Kingman's theorem, with probability one, there exists T^(x) such that if t > T3, then 

d(0,tx) > (1 - 6)tn(x). 

Choose finitely many Xi G 2A \ (1 + e)A such that the closure of 2A \ (1 + e)A is covered by (J B E (xi, 5e 2 ). 
Belonging to 2A \ (1 + e)^4 implies that ^(ajf) > (1 + e). Let T = max i {T 2 (a;i), T 3 (xi)}, and let t>T. Then 
ifyeB^Si.ae 2 ), 

d(0,ty) > d(0,txi) — d(txi,ty) 

> (1 - 5)tfj,(Xi) - e 2 t 

> {1 - S)(l + e)t - e 2 t 

> {1 + e 2 )t - e 2 t = t, 

where the final inequality follows from the assumption (|3.3[) on <5. Thus (J B E (xi, Se 2 ) C i-BJ 1 so 

2A \ (1 + e)A C |J B E (x l; 5e 2 ) C lB t c . 



The set jB t is connected and contains the origin, hence jB t C (1 + e)^4 as desired. 



□ 



Corollary 3.2. If the measure P is isotropic (rotationally-invariant), then the limiting shape A is the 
Euclidean ball of radius □ 

A simple consequence of Lemma l2.4l is that the convergence (|1.3[) given by Kingman's theorem is uniform: 

Proposition 3.3. For all e > 0, with probability one, there exists T > such that if t > T, then for all 

I jd(0, tv) — fi v \ < e. 

Proof. Suppose not. Then with positive probability, there exist e > 0, t n — > 00 and v n G S' ci_1 such that 



j-d(0,t n v n ) - jj, Vfl 



> e, 



for all n. By compactness of the sphere, a subsequence of v n converges to some v G S 1 ; assume without 
loss of generality that v n —¥ v. By the triangle inequality, 



(3.4) 



e < 



±d(Q,t n v n ) - ±d(0,t n v) + ±d(0,t n v) - fi v + \fj, Vn - fj, v \ 



The second and third terms tend to zero a.s. as n — > 00 by Kingman's theorem and Proposition 11.31 
respectively. By the triangle inequality, the first term is bounded by j-d(t n v,t n v n ). Let p = e/2K, where 
K is as in Lemma 12.41 Since for large n, \v — v n \ < p, we can apply that lemma with x = v to get 



T^d(t n v, t n v n ) < Kp = e/2, 
for large n almost surely. This contradicts 



□ 



Theorem 3.4. With probability one, if 7 is a C 1 -curve parametrized by Riemannian length, then \ j(t)\ < 00 
for all t > 0. 
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Proof. It suffices to consider curves starting from the origin. Suppose that with positive probability, there 
exists a smooth curve 7 parametrized by Riemannian length which starts from the origin and for which 

(3.5) lirn| 7 (i)|=oo 

for some finite T. Since 7 is parametrized by Riemannian length, for alH > 0, 

t>d(0, 7 (t)). 

Proposition [33] and (13.5[) imply that 

„ ,• t ,• d(0,l(i)) 

= hm - — r— > hm — ; — — — > mmu^ > 0, 

t^T \j(t)\ ~ t->T | 7 (i)| " v 

a contradiction. □ 

We can refine this result if we impose an additional smoothness constraint on the metric. If g is a C 2 - 
smooth Riemannian metric, i.e. g <E C 2 (R d ,SPD), then we can use the calculus of variations to derive 
the Euler-Lagrange equations for the functional R. These are called the geodesic equations [13] for the 
Riemannian metric g, and solutions to this system are called geodesies. The geodesic equations form a 
second-order system with locally-Lipschitz coefficients, so a geodesic is uniquely determined by its starting 
point and velocity. We call a geodesic 7 length-minimizing if for all x, y € 7, the distance d(x, y) is realized 
as the Riemannian length of the part of the curve 7 which connects the two points. Not all geodesies are 
length-minimizing; for example, on the sphere, the geodesies are great circles, which do not minimize length 
past antipodal points. 

A metric is said to be geodesically complete if for all x £ M. d and v £ T x M. d , the unique geodesic 7 at x in 
direction v can be continued for all time. Part of the Hopf-Rinow theorem |13j of Riemannian geometry is 
that geodesic completeness is equivalent to the condition stated in Theorem 13 .41 for our random metric g. A 
further corollary |13) is that distances are always realized by geodesies. Summarizing, we have 

Corollary 3.5. Suppose that, in addition to the assumptions of this paper, g is a C 2 -smooth random 
Riemannian metric. Then, with probability one, g is geodesically complete. Consequently, for all M. d , 
there is a finite length-minimizing geodesic 7 connecting x to y such that 

d{x,y) = R{ 1 ). 

We alert the reader to a different meaning of the word "geodesic," used often in the first-passage percolation 
literature. The term is used there to denote a globally length-minimizing path. This is very different from 
the standard meaning of the word in differential geometry: as described above, geodesies arc the curves 
which locally minimize length, but not necessarily globally. We adhere to this meaning in the present paper. 

The existence of two-sided minimizing paths is an open question for all FPP models. For two-dimensional 
standard FPP, Licea, Newman and Piza [TH [TSl US] have a number of results in this direction. Their work 
relies on certain curvature assumptions about the limiting shape, which have not been verified for models of 
independent FPP. Howard and Newman's j9][T0] model of Euclidean FPP, on the other hand, is rotationally 
invariant. Consequently, the limiting shape is a Euclidean ball, and they prove many results not available 
in the lattice setting. See the excellent survey [8] for more details. By Corollary 13.21 the limiting shape for 
isotropic Riemannian FPP is a Euclidean ball. Our hope is for this setting to be a fertile ground for adapting 
the results referenced in this paragraph. 

Remark added after publication: The moment assumptions (|1.2[) and (|2.1[) can be slightly weakened 

to: 

M (r) = E[e rA °] < 00 for all r < a, 
for some a > 0. The only modification to the paper is in the proof of Lemma [2~2l where one uses M — M(a) 
instead of M = M (1). We leave the details to the reader. 
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